arXiV: hep-ph/011215 2vl 11 Dec 2001 


UT-978 


Proton Decay in the Semi-Simple Unification 


Masaaki Fujii and T. Watari 

Department of Physics, University of Tokyo, 
Tokyo 113-0033, Japan 


Abstract 

Semi-simple unification is one of a model which naturally solves two difficulties 
in the supersymmetric grand unification theory; doublet-triplet splitting problem and 
suppression of dimension 5 proton decay. We analyzed the dimension 6 proton decay 
of this model using perturbative analysis at the next-to-leading order. The life time of 
proton is 3 X 10^^ - 10^^ years for wide range of SUSY breaking parameters, and there is 
an intriguing possibility of observing proton decay signals in the next-generation water 
Cerenkov detectors such as Hyper-Kamiokande and TITAND. Several uncertainties in 
this prediction are also discussed. 


Introduction 

Supersymmetric(SUSY) SU(5)gut grand unification theories(GUT’s) are supported by 
the approximate SU(5)gut unification at around 10^®GeV of the three gauge coupling con¬ 
stants of the minimal supersymmetric standard model(MSSM)[2]. However, the conceptual 
beauties of the GUT’s[0 and such a phenomenological success are not more than indirect 
evidences, and it would be the proton decay signals that makes us believe the GUT’s in 
nature. 

The minimal SU(5)gut rnodel predicts proton decay through dimension 5 operators^, 
and is now almost excluded^ by experimental bounds such as r(p —^ K^v) ^ 6.7 x 10^^yr.(90 
% G.L.)|^. Therefore, we have to analyze the proton decay in an extended model in which 
those operators are suppressed or absent. Predictions on the proton decay through dimension 
6 operators severely depend on how a model is extended^; the life time of the proton depends 
on the fourth power of the mass of SU(5)GUT-off-diagonal gauge boson (GUT gauge boson), 
and hence on the detailed spectrum of the model around the GUT scale. Therefore, an 
analysis on the proton decay has to be based on a phenomenologically reliable model of the 
GUT’S. 

The doublet-triplet splitting problem 0 and suppression of the dimension 5 proton decay 
operators had been the two major obstacles in model buildings of the GUT’s. The semi¬ 
simple unihcationfl, P, is a model that solves these two problems in a natural way. In this 
letter, we calculate the proton decay rate in this model. The proton decay is relatively fast in 
this model, whose reason will be clear in the text. We restrict ourselves in parameter region 
of the model where a perturbative analysis is valid. As a result of a full next-to-leading order 
analysis[0], we found that the mass of the GUT gauge boson can be determined, and that 
the resulting life time of proton does not depend on SUSY breaking parameters so much: 
the life time is r ~ (3 — 10) x 10^'^yr. Various sources of uncertainties in this prediction are 
summarized at the end of this letter. This result means that the proton decay is generically 
detectable in the next-generation water Cerenkov detectors such as Hyper-Kamiokande and 
TITAND 0, g. 

Brief Review of the Semi-Simple Unification Model 

We briefly review the semi-simple unihcation model that uses SU(5)gut x U(3)h gauge 
group. This gauge group is directly broken down to the SU(3)c x SU(2)l x U(1)y of the 
MSSM. Quark and lepton (5* -f 10) and Higgs (i7*(5) + Hi{^*)) supermultiplets are singlets 
of the U(3)h gauge group and transform under the SU(5)gut as in the standard SU(5)gut 

^ For example, in some type of model , the dimension 6 operators are not induced by the GUT gauge 
boson exchange. 
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model. Some more fields are required for the GUT symmetry breaking: X“g(Q;,/3 = 1,2,3) 
transforming as (l,adj.=8 + 1) under the SU(5)gut x U(3)h gauge group, and = 
1, • • •, 5)+Q^q and Q\{i = 1, • • •, 5)+Q^^ transforming as (5* + 1 , 3 ) and (5 + 1 , 3 *). Indices a 
and f3 are for the U(3)h and i for the SU(5)gut- is also expressed as X'^(fc)“/ 3 (c = 0, 1 — 8 ), 
where t°‘{a = 1 — 8 ) are Gell’mann matrices^and to = Isxs/n/S- Superpotential is given asHTOj 

-ylAmwW^ ( 1 ) 

+h'm\Q% + 

Tj/iolO ■ 10 ■ H + 2/5*5* ■ 10 ■!? + ■••, 


where the parameter v is of order of the GUT scale, ywiV^* are Yukawa coupling constants 
of the quarks and leptons, and Ash, Ag^, Am, Agjj, h', h are dimensionless coupling constants. 
One can see that the above superpotential has Z 4 R symmetry under a charge assignment 
given in Table |I|, and this symmetry forbids enormous mass term W = HH for the Higgs 
doublet^. The bifundamental representation Q" and Q\ acquire vacuum expectation value, 
(Q“i) = and { Q\ ) = because of the second and the third line in ([|). Then, 
the mass terms of the colored Higgs multiplets arise from the fourth line in (|I]) in the GUT- 
breaking vacuum. The Z 4 R symmetry is not broken even after the GUT symmetry is broken. 
One can also see that this Z 4 R symmetry forbids the dangerous dimension 5 proton decay 
operators lU = 10 ■ 10 ■ 10 ■ 5*. 

Naive Estimation 

At tree level, the gauge coupling constants of the SU(3)c x SU(2)l x U(1)y are given in 
terms of those of the SU(5)gut x U(3)h as : 


03 acJ ogut ctsH’ 


( 2 ) 


and 



CX\ (Xy ) 


1 

ogut 
1 2/5 

- 1 - 5 

OgUT CHH 


( 3 ) 

( 4 ) 


normalization condition tr{tatb) = (5ab/2 is understood. Note that the normalization of the following 
to is determined so that it also satisfies trfioA) = 1/2- 

^/i-term can be obtained through the Giudice-Masiero mechanism [p^. 
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where ac, wl, «gut, «3H and am are fine structure constants of the three MSSM gauge 
groups, SU(5)gut, SU(3)h and U(1)h0, respectively. The approximate SU(5)gut relation and 
deviation form it (Fig.^) are naturally explained through the above equations if I/ogut ~ 24 
and l/a 3 H ~ 1, l/aiH ~ 2.5. At the same time, we notice that the “GUT scale”, an energy 
scale at which Eqs.(|]-§) hold, lies lower than the ai - 02 unification scale Mi _2 and higher 
than the 02 - a^ unification scale M 2 - 3 . Therefore, the decay rate of proton is expected to 
be enhanced compared with the rate using the Mi _2 as the GUT gauge boson mass. 

At one-loop order, the gauge coupling of the U(1)h runs0 asymptotic non-free : 

There are two important remarks here. First, the cut-off scale A of this model exists below 
the Planck scale; 1/am should be positive even at the A. The constraint 1/am^ 2.5 at the 
“GUT scale” allows the A to be higher than the “GUT scale” by one order of magnitude or 
a little bit more, and this is much enough to justify the field theoretical description of the 
GUT symmetry breaking and to accommodate all the GUT spectrum below the cut-off scale 
A. This A lies around lO^^GeV or a little higher, though it may be below 10^®GeV. Secondly, 
the IR-free (asymptotic non-free) behavior of the U(1)h coupling leads to 


1 1 

-<- 

ttSH ttlH 

at the “GUT scale” under an assumption 


( 6 ) 


—(A) ^ —(A), (7) 

ttSH (y-lH 

This U(3)h- relation at A is quite natural if there is U(3)H-structure in a fundamental the¬ 
ory [^. Then, as a consequence of the relation Eq.(|l), we notice that the “GUT scale” is 
closer to the M 2-3 rather than to the Mi _2 since (1/ai — 1 /q! 2 ) > (l/as — l/a 2 ) there, and 
hence the proton decay is relatively fast. 

Threshold Corrections at the GUT Scale 

In the analysis at the next-to-leading order, one-loop threshold corrections of the GUT 
model are also taken into account. The three MSSM gauge coupling constants just below 
the GUT scale are expressed in terms of the gauge coupling constants and various masses in 
the spectrum of the GUT model, including the mass of the GUT gauge boson, M^. Particle 

^The normalization of the U(l)H-generator is determined so that Q, Q have charge ±1/V6. 

^The one-loop /3-function of the SU(3 )h coupling is 0. 
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spectrum around the GUT scale is summarized in Table Explicit expressions of the MSSM 
gauge coupling constants are given as follows : 


1 

as 


(h) 


1 

02 


(h) 


-(A) + 

OgUT 


^(A) + ^ln 
asH ^ ^ 27r 






OgUT 



— In (—) 
271 \Myj 


6 , 

H-In 

271 


Msy\ 

Wc) 






( 8 ) 

(9) 


1 


(h) 


-(A) + 

OgUT 



am 



277 



2 jr \MiMJ 




( 10 ) 


where /x is a renormalization point, which is taken to be just below the GUT scale, My, My, 
My, Mgy, Mgc are masses of particles around the GUT scale (see Table |^) and aGUT, 3 H,iH(A) 
are hne structure constants of the gauge groups SU(5)gut, SU(3)h and U(1)h, respectively, 
at the cut-off scale A. 

In general, it is impossible to determine the My if GUT models have more than three 
parameters. However, it is not necessarily the case in the semi-simple unihcation model. 
Threshold corrections in Eqs.(|8|-pTD is simplihed considerably under two reasonable assump¬ 
tions. One is the U(3)H-relation Eq.(0) and the other is Af = 2-SUSY-relation: 


9m — Aih(~ Ajjj), 


5'3H — A3h(~ Agy). 


( 11 ) 


Under the latter condition, a large threshold correction from the massive SU(3 )h vector 
multiplet^ is almost canceled by its Af = 2-partner, the SU(3)c-adj. chiral multiplets, since 
Mgy ~ Mgy. Now that the threshold corrections form the SU(3)c-adj. multiplets decouple 
from Eqs.f^pUl), we are left only with two threshold corrections: one from the massive vector 
multiplet of the GUT gauge boson and the other from colored Higgs chiral multiplets. Then, 
one can easily see that three combinations. 


1 

aGUT 


‘“(4 


and 


1 

asH 


(A) + -l„ 



( 12 ) 


are determined in terms of the values to be put in the LHS’s and deviation from the U(3)h- 
relation and the Af = 2-relation, once the cut-off scale A is hxed. 

In particular, the GUT gauge boson mass is given by 



®Note that Msv ~ 





10 • My, and hence the threshold correction is large. 




(13) 
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The last two factors show how the result is changed due to the deviation from the assumptions 
we made. A“^/^-dependence is a direct consequence of the one-loop running of the am in 
Eq.(|^), and this negative power dependence implies that this gauge boson mass is generically 
light. The life time of proton through this GUT gauge boson exchange is given in terms of 


the My as|ll 


t{p —>■ TT^e"'') ~ 0.61 X 10 


35 


X 


(- 


M, 




VlOieCeVy l24aGUT(My) 


'0.15(GeV)2 


yr., 


(14) 


where lU is a hadron matrix element calculated with lattice quenched QGD||l6|. 

Threshold Corrections at the Weak Scale and Two-Loop Rnnning 

In order to determine the precise value of the GUT gauge boson mass by using (|^), 
we must accurately determine the three MSSM gauge coupling constants at the GUT scale. 
For this purpose, we take full one-loop threshold corrections at the weak scale into account 
for the three gauge coupling constants and top- and bottom-Yukawa coupling constants by 
following the method in Ref. and use two-loop renormalization group(RG) equations. 
For illustration, let us briefly review the procedures which we adopt in this letter. The 
conventions of SUSY breaking parameters and of the sign of the p-term are the same as 


those in Ref. [17 


The SUSY threshold corrections to £x the DR coupling a^{Mz) is very simple and the 
result is as follows: 

“3(^z)mS 


o^siMz) — 


1 — Aaa 


where 


Aa-j = 




27r 


-In 


mt 

2 


21n 


m-g 

Mz 


EEln 


i=l 


m, 


Qi 


Mp 


(15) 


(16) 


Here, Mz = 91.188GeV is the Z-boson pole mass and we take as^Mz)^;^ = 0.118(2) p0| . 
The summation with q runs over all the six squark flavors, and the constant contribution 
1/2 in Eq. (p[6D is necessary when the coupling is translated from the MS scheme to the DR 
scheme. 

Because of the breaking effects of the SU(2 )l gauge group, the determinations of the 
DR gauge coupling constants ay^Mz) and ai^{Mz) are much more complicated. First, we 
calculate the DR electromagnetic coupling constant, a{Mz)- The explicit formula is given 
by 


Oi{Mz) = 




1 


1-Aa’ 


Oif^ra — 


137.036’ 


(17) 


5 




























where 


Aa = 0.0682 ±0.0007 
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f mH+ \ 4 


^ d i=l 


\Mo 


+ oEEin 
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+ 9 EE 1 “ 
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, 4 A, f'^X+ 
+ Eln' ^ 

i=i 


mu.^ 

M, 


My 


.(18) 


Here, J2u denotes a sum over u, c, t, and similarly for J2d S'lid J2e- The numerical values ap¬ 


pearing in the above expression includes the two-loop QED and QCD corrections in Ref. []^ 
as well as the hve-flavor contributions in Ref. 0- 

Next, we need to £x the DR weak mixing angle 6*ew to derive the DR gauge coupling 
constants, ayi^Mz) and a\^{Mz)- The formula to get the DR weak mixing angle is given by 

na^Mz) 


cosR6'ew)sinR6'ew) = 


Ar = 


V2M|G^(1-Ar) ’ 


P 


n^ty(O) 


TZe 


nE(M: 


±h 


VB 


(19) 


( 20 ) 


where Mw = 80.419GeV is the W-boson pole mass, p is defined as p = M^/(cos^(6'ew)Tf|), 
= 1.16639 X 10“®GeV“^ is the Fermi constant, and Sys denotes the nonuniversal vertex 
and box diagram corrections. The explicit formulae to calculate the quantities given in the 


above expressions and the DR Yukawa coupling constants are all given in Ref. [0 


Taking ttem, ol^{.Mz)^.i the quark and lepton masses, and SUSY particle masses as inputs, 
we calculate all the three gauge coupling constants and top- and bottom-Yukawa coupling 
constants in DR scheme at the Z-boson pole mass with full one-loop threshold corrections. 
With these values and tree level tau-Yukawa coupling, we use the two-loop RG equations to 
obtain the gauge coupling constants at the GUT scale. For the Yukawa coupling constants, 
we use one-loop RG equations. 

In this letter, we adopt the central values given in Ref. |2^ for the masses of vector bosons, 
quarks and leptons[|. As for the mass spectrum of the SUSY particles and light Higgs particle, 
we take the values calculated by the SOFTSUSY code 0| with mSUGRA boundary condi¬ 
tions for demonstration.^ By using these input values with mSUGRA boundary conditions, 
we also conhrm that the unihcation-scale correction 


<y3{Mi-2) — ai(Mi_2)(l ± 2eg), 


( 21 ) 


^Neutrino masses are set to be zero. 

® We greatly thank K.Suzuki for this task. 
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at the ai -«2 unification scale Mi _2 is quite consistent with the result given in Ref. |]^ . 

Conclusion 

Now, we can estimate the proton life time for various SUSY particle spectra. We neglect, 
for the moment, possible two uncertainties expressed by the last two factors in (|T^) coming 
from the deviation from M = 2-relation and U(3)h- relation. Effects of such deviations are 
discussed later. Here, we also set the cut-off scale to be lO^^GeV; In most part of SUSY 
breaking parameter space, the three gauge coupling constants unify approximately at around 
lO^^GeV and hence the cut-off scale A is expected to be no less than lO^^GeV. Therefore, 
we obtain a conservative upper bound of the proton life time, using the lowest cut-off scale 
(see (0)). 

We plot the contours of the life time of proton in mo — M 1/2 plane, where mo and M 1/2 
are the universal soft scalar mass and gaugino mass at the GUT scale, respectively. In 
Fig.§, we show contours of the proton life time for /i < 0 cases with several choices of 
^0 (= 0, —300GeV), the universal A-term at the GUT scale, and tan/? (= 10, 30). The 
contour plots for /i > 0 cases are given in Fig. |^. 


As we can see from these contours, the proton life time is in the range 3 x 10^“^ — 10^® yr. 
in most part of the parameter space regardless of choices of tan/?, Aq and sign of p. We find 
the minimum of the proton life time is no less than 3 x 10^“^ yr. in whole parameter space, 
which is well above the current experimental limit by the Super-Kamiokande, 5.0 x 10^^ yr. 


(90% G.L.)|^, 1^. The thick gray contour lines corresponding to the life time of proton 
7 X 10^“^ yr. represent the 3 ct discovery limit of the IMt (fiducial volume) detector after ten 
years running |U 


Therefore, in the semi-simple unihcation model, we have an intriguing possibility to con- 
hrm the existence of the GUT in nature by observing the proton decay in the next-generation 

In the op- 


Mt water Gerenkov detectors, such as Hyper-Kamiokande and TITAND 
timistic cases with some enhancement factors of the decay rate of proton (see below), we have 
a chance to detect the proton decay also in UNO (~ 500kt hducial volume) experiment. 

Although we set the cut-off scale A to be lO^^GeV in calculating the GUT gauge boson 
mass to obtain the conservative lower bound of the proton decay rate, the actual cut-off 
scale may be a little more higher. In that case, the rate is enhanced by (A/lO^^GeV)^. 
Another possible enhancement of the decay rate arises when there are SU(5)GUT-charged 
particles at an intermediate scale. Existence of such particles are highly motivated in the 
semi-simple unihcation model; 5-1-5* representations are required at the TeV scale when the 
discrete Z 4 R symmetry is gauged since the discrete gauge anomaly Z 4 r-[SU( 5 )gut]^ should 
be canceled [^. In this case, the gauge coupling constant ogut is stronger as a result of 
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the RG flow with new particles, and the decay rate is enhanced by 1.6. Although one might 
suspect that there is a one-loop threshold correction from a possible mass splitting between 
triplets and doublets in 5 4-5*, and that the GUT gauge boson mass would be also changed, 
the GUT gauge boson mass is actually stable against this correction, since Eq.(|^ is an 
expression from which the threshold corrections from the colored Higgs multiplets decouple. 
The same thing happens when the SUSY breaking is mediated through gauge mediation 
because of the presence of the messenger sector, though the SUSY threshold correction should 
be re-analyzed using the spectrum of the gauge mediated SUSY breaking in that case. 

Finally, we summarize various uncertainties in the theoretical prediction given above. 
The first uncertainty comes from possible violation of the U(3)h relation. The violation 
|(1 /q!3h — l/aiH)(A)| = 1/3 leads to a change in the decay rate by x/ 4- 0.5. The second 
uncertainty comes from an error bar of the experimental values of the QGD coupling. This 
results in uncertainties by factor x/ 4- 0.7 for la error. The calculation of hadron matrix 
element in has an error W = —0.153(19)GeV^, which leads to a factor x / 4-0.8. Another 
uncertainty comes from a possible non-renormalization operators involving the {QQ^ vacuum 
expectation value in the gauge kinetic function of the SU(5)gut0 • They generically modifies 
the SU(5)gut relation directly by {QQ^ /A^ ~ 10“^ at tree level. If it is the case, the possible 
change in the result will be at most roughly the same amount as those discussed above. 

There are two more sources of uncertainties whose effects we cannot estimate. First, if one 
considers an exotic situation in which unknown non-renormalizable operators are relevant in 
the Wilsonian RG equations, then the perturbative analysis we adopted in this letter is not 
adequate since we omitted such effects. Secondly, we cannot estimate anything without the 
M = 2-SUSY relation. This is because the perturbative analysis above the GUT scale is no 
longer valid without this relation, as is discussed in the appendix. 
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® Such non-renormalizable terms in the gauge kinetic function is expected to be suppressed when one 
considers a certain structure of the fundamental theoryfl^. 





A Role of approximate J\f = 2 SUSY relation 
in perturbative analysis 


The GUT-breaking sector of the semi-simple unification model has a multiplet structure of 
M = 2 SUSY, and the interactions between them (the hrst - the third lines in Eq.(|I])) are 
quite similar to the M = 2 gauge interactions with Fayet-Iliopoulos F-term. Therefore, it 
is quite likely that this apparent M = 2 structure is a remnant of the M = 2 SUSY in a 
fundamental theory|T5[|. Then, the approximate M = 2 relation Eq. (|TT|) at the cut-off scale 
would be a natural consequence. 

The approximate M = 2-relation is not only expected as above, but also almost required 
from another reason. The perturbation analysis performed in the text is no longer valid if it 
is not satished and that is the reason why we assumed this relation throughout this paper. 

Let us suppose that the couplings Ag^ and in the superpotential (|I]) are large compared 
with ^fgH and ^^ih. Then, those couplings become large extremely fast through one-loop RG 
equations, and hence we have to require that agjj = Agjj/(dvr) and = A^jg/(^^r) are well 
below 2a3H and 2aiH, respectively. The same discussion also holds for Ag^ and A'^y. Now 
what if those couplings are small compared with the corresponding gauge couplings? In this 
case, we can neglect the last two terms in the following two-loop RG equations of the gauge 
couplings. 


d 

clln/i 

d 

clln/i 



Qih + 17q;3H |(q^ih + , |(«ih + 

27r2 27r^ 2 ^“^ 

_6_ Qih + 8q3h |(q^ih + Sasn) |(q^ih + 8«3h) 

27r 27r2 2 ^“^ 2 tt ‘^ 


( 22 ) 

(23) 


Then, am becomes large quite rapidly and am becomes large more faster than in the one- 

(^^) 

loop running. Thus, we require that and comparable to the gauge couplings so 

that the two-loop effects are negligible. 

In the approximate M = 2-SUSY limit and only in this limit, am — and am — 'aiH\ 
anomalous dimensions of hyper multiplets. 


iQi 


+ <aiH _ SftgH + am 
Gtt Gtt 


(24) 


vanish at all order, and the RG flows of the gauge couplings are one-loop exact. Then, in 
turn, all other parameters in the superpotential, in particular h and h', are stable against 
quantum corrections from the strong couplings am, « 3 H, and 
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Values of the coupling constants h and h' themselves are the possible obstruction left 
behind for the perturbative analysi^. They are obtained from a ratio y/M^Mc/M^, which 
in turn is obtained from Eqs.(|8|-pTD in the way described in the text: 


Vhh' = v^^fQUTe “3H 




(25) 


The value of the RHS of this equation varies from sub-(!2(l) to 0{1). Therefore, we can 
expect that the perturbative analysis performed in the text is valid for most part of the 
SUSY breaking parameter space, taking into account the uncertainties in the gauge coupling 
constants. 
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Fields 

5*,10,1 

H,H X-g Q,,g* 

Qe 

Z 4 R charge 

1 

0 2 0 

2 -2 


Table 1: Charge assignment of the Z 4 R-symmetry is given. 1 denotes a right handed 
neutrino. 


(3,2)-t 



(i,ir 

(1,1)° 

(adj., 1)0 

(adj., 1)0 

m.vect. 

x + x"^ 

X + X^ 

m.vect. 

X + X^ 

m.vect. 

X + X^ 

M,= 

M,= 

M,= 

Miv = 

Mi,= 

Msv = 

Msc = 

V^gGVTV 

hv 

h'v 

+ 25'gut/5)'^ 

V^AihU 

\/2(5'|h + 5'gut)'^ 

y/2X'^HV 


Table 2: Summary of the particle spectrum around the GUT scale. The hrst line denotes 
the representation under the MSSM gauge group. In the second line, m.vect. denotes Af = 1 
massive vector multiplet and X + a pair of A/” = 1 chiral and anti-chiral multiplet. Mass 
of each multiplet is given in terms of gauge couplings and parameters in the superpotential 
(|l|) in the fourth line, and given in the third line is the expression of the mass used in the 
text. Multiplets with masses Mi^ and Mi^, and Mgc can be regarded as J\f = 2-SUSY 
partner with each other in the Af = 2-SUSY limit (see also appendix). 
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Figure 1: Approximate SU(5)gut relation between the three MSSM gauge coupling constants 
and deviation from it. Icr error bar of the QCD coupling are also described. SUSY threshold 
corrections are calculated using the spectrum of mSUGRA model with tuq = 250GeV, M 1/2 = 
SOOGeV, Aq = 0 and tan/3 = 10. The sign of /i-term is taken to be negative. 
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tan/3 = 10, ^0 = 0 tan/3 = 10, ^40 = —300GeV 




tan/3 = 30, = 0 tan/3 = 30, ^40 = —300GeV 




Figure 2: Contour plots of the proton life time in mo — M 1/2 plane for /i < 0 cases. SUSY 
threshold corrections are calculated using mSUGRA sparticle spectrum with universal bound¬ 
ary conditions mo, M1/2, Ao at the GUT scale. As for (tan/3, Ao[GeV]), we take them to 
be (10,0), (10, —300), (30,0), (30, —300), respectively as you can see from each hgure. Solid 
lines correspond to the contours of the proton life time, 5 x lO^V-, 7 X lO^V-, lO^V-, 
2 X lO^^yr. from in to out, respectively. Some of them are explicitly denoted in each hgure. 
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tan/3 = 10, ^0 = 0 tan/3 = 30, Aq = 300GeV 




Figure 3: Contour plots of the proton life time in mo — M1/2 plane for /i > 0 cases. Other 
conventions are the same as those in Fig. |^. 
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